Thermoelectric efficiency of heterogeneous media at low temperatures 



o 
o 



> 



1^ 
o 



X 



S.A. Ktitorov, V.K. Zaitsev, M.I. Fedorov 
A.F. loffe Physico-Technical Institute of the Russian Academy of Sciences, 
Polytechnicheskaja str. 26, 194021 St. Petersburg, Russia 

' A quantum limit of the thermoelectric efficiency for heterogeneous media with weak links is 

I established with a use of the Landauer-type formulae. 
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^ \ I. INTRODUCTION 

^: . 

. As is well known, the kinetic parameters of semiconductors and metals such as conductivity a, thermoelectric 

' power S and thermoconductivity k are determined by a crystallographic and atomic structure of the material and 

] are, therefore, very specific for an every given material [l|. However, there is a class of conductors that have low- 
temperature kinetic properties, which are not specifically related to a concrete scattering mechanism and to other 

Q , nonuniversal factors. We will show here that the thermoelectric efficiency zT — ^ of these structures at low enough 

D ■ temperatures can be expressed in terms of the fundamental constants and of the system geometry. 

a: 

^ • II. CONDUCTING STRUCTURES WITH WEAK LINKS 

> i 

I ; ' Let us consider a polycrystal with point contacts between crystallytes. These point contatcs are really small 
d bridges or weak links between the crystallites. We assume that within the actual temperature range, all quasiparticles 
(charge carriers, phonons, etc) are almost thermalized in crystallites and their distribution functions insignificantly 
I , differ from the equilibrium ones. Large volume of crystallites in comparison with bridges connecting them ensures 
' ^ ■ effective thermalization. Notice, that by our opinion, it is not even necessary for strong inequality D 3> I (with I 
being the mean free path) to be fulfilled because nonspherical crystallites can play a role of three-dimensional Sinai's 
billiards '^l for most of the quasiparticles. In other words, large time delay in propagation of quasiparticles in result 
of multiple reflections on the crystallite boundary leads to effective thermalization, if the total path through the 
crystallite including multiple reflections by the crystallite boundaries exceeds the mean free path. Similar processes 
can occur in conducting composites based on opals, porous glasses etc. Rather narrow windows in the spectrum are 
formed by excitations, which avoided multiple reflections because of a spatial coherence resonance. Such windows 
ly-^ of coherence can be actual only at very low temperatures. Large cross section of a crystallite makes the electric 
\^ , current and heat flux densities extremely small there. Thus, we can consider the system as a network of electric 

■ and thermal resistors and thermoelements (weak links) connecting nodes (crystallites), which can be characterized 
' by their electric potentials and temperatures .Then the transport problem is reduced to the following two steps: (i) 

calculation of the electric and thermal resistance and of the thermoelectric power of a single weak link; (ii) calculation 
of integral characteristics of the network. The latter is extremely simple in the case of a regular network, but can be 

■ rather complicated in the irregular case. In the case of a strong irregularity (severed connections, in particular), the 
' percolation theory can be applied. 



^ • III. TRANSPORT PROPERTIES OF A WEAK LINK 

o : 

I Let us consider a bridge connecting two crystallites as a figure of revolution only slightly different from a cylinder. 
^ ' Phonon and electronic contributions have to be considered separately. 

Long wavelength excitations propagate without feeling of possible disorder that allows us to use the quasiparticle 
picture at low temperatures. Size quantization of excitation spectra in the weak link permits us to consider it as a 
' waveguide for excitations that establishes a lower limit for propagating waves using the relation: 



/?„ - ^/k^~Tcf„ (1) 

where k is the wave vector value, /3„ is the wave vector component along the waveguide axis, x„ is the membrane 
eigenvalue, the lowest one is of order d~^. Simple estimates show that these wavelengths are related for acoustic 
phonons and electrons, respectively, to energies of order 1 K and 10 K. Therefore, only a few waveguide modes 
participate in transport processes within this temperature range, but their number fast increases with energy. This 
conclusion is valid if we assume that actual excitations feel the weak link as a regular waveguide that can be correct 
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only for the most low-lying excitations with large wave length. When the energy of an excitation increases, the bridge 
becomes more irregular for waves. The transmission coefficient decreases catastrophically except for a few (if they 
exist) special energy values related to possible resonances in irregular interfaces. Besides, there exist gapless surface 
modes. 

The one-dimensional electronic ballistic transport in the presence of the temperature and electrochemical potential 
is described by the equation 



2e^ 



E 



dE 



dnpiE) 
dE 



t:\E) 



E- fi 
T 



AT + A^l 



(2) 



where J is the electric current through the link, np is the Fermi-Dirac distribution, tfj' (E) is the electron transmission 
coefficient, fi and A/i are respectively the electrochemical potential and its drop across the weak link, AT is the 
temperature difference between the adjacent grains. The energy flow Q = Q^^ -I- Q"^' comprises contributions both 
from the electronic and phonon subsystems 



where (fc) is the v-Vn acoustic phonon branch frequency, n^'' and rJ'l^ are respectively the quasiequilibrium phonon 
distribution functions in the bulk of the right and left crystallites; e^^ [k] is the v-th branch of the electronic spectrum 
in the link, .n^ and nf^ are corresponding electronic distribution functions. 

IV. KINETIC PROPERTIES OF THE NETWORK 

Within the approximation of locally equilibrium crystallites, the electric potential and temperature distributions 
are governed by a finite-difference equations. In the steady state, the electric potential and temperature distributions 
satisfy the time- independent finite-difference electric and heat transport equations: the "Kirhgoff law" for currents 
entering into the n-th site 

= 0, (5) 

(m) 

^Qnm = (6) 

(m) 

and the Ohm-Fourier-Peltier-Zeebeck law for the electric Jmn and heat Qnm currents between the sites m and n 



where n stands for the crystallite number vector with components n^, Uy, n^; Gnm is the thermal conductance of 
the weak link connecting crystalites m and n, Gmn 7^ if and only if m and n differ by a unit superlattice vector 
e.; (m) means summing over the nearest neighbours of the given node n., i?nm is the electric resistance of the link, 
Vn is the potential of the n-th crystallite, the thermopower Snm and the Peltier coefficient Ilnm are related by the 
Kelvin-Onsager symmetry law Ilnm = SnmiT^ + ?m)/2, Qnm is the energy flow through the link, Jnm is the electric 
current through the link. Now the theoretical analysis of the heat transport problem for the polycrystal is reduced to 
two different ones: the first problem is to study the electric and heat transport through the weak link and the second 
is to solve the network equations for a given distribution of links parameters. The former can be considered using 
different methods depending on the temperature range, the weak link geometry, and the microscopic properties of the 
material, while the latter is determined exclusively by the system geometry and a distribution of parameters. For a 
regular superlattice we have Gn.n+e = G for all n and m, and the effective thermoconductivity of the sample reads 



>^efS = G/A 



(8) 



where A is the superlattice spatial period. A calculation of the effective thermoconductivity in the case of a random 
distribution of Gmn can be very sophisticated problem mathematically similar to the electron transport in irregular 
conductors, but we concentrate our attention here on the weak link thermal resistance assuming that the superlattice 
is regular. 



V. EFFICIENCY 



Assuming the structure to be regular, we calculate the thermoelectric efficiency of a single link. The thermoelectric 
efficiency can be expressed in terms of the single link transport coefficients: 
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where S, Gei, and Gth are respectively the thermopower, the electric and thermal conductances of the single link, 
which are given by the formulae following from [21 ; 



fee ■\p 
e ^ 



dE ■ t 



,1 dn"^ E 



" dE kuT 



'^<'iE 



(10) 



2e^ 

2'Kh 



(11) 



E 



3xexp(a:) 

2 V {■''^J 



dx 

[exp (x) - 1] 



(12) 



Notice that the group velocity ^ and the one-dimensional density of states 1/ (f|) cancele one another in the 
one-dimensional transport. Thermal conductance is extremely sensitive to the real geometry of the link that can 
be seen from the effective dimension of the transport and the behaviour of the transmission coefficient. In the 
three-dimensional diffraction regime it reads: 



(13) 



In the one-dimensional case it can be written as 



■TN, 



(14) 



where a is the link cross-section area, N^^ is the number of the gapless phonon modes. 
We can write within this approximation for the electric conductance : 



G"' = ^N^\ (15) 

While the conductance and the thermoconductance have a universal form, which is very not sensitive to details of the 
weak link geometry, 
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Substituting these formulae into[9l we obtain for the thermoelectric efficiency quantum limit at the low-temperature 
regime 



This means that one has to minimize the number of modes involved into the transport through the weak link in 
order to rise the thermoelectric efficiency. 
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